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We study the coexistence of strain- and charge-mediated magnetoelectric coupling in a cobalt 
(0—7 nm) wedge on ferroelectric [Pb(Mgi/3/Nb2/3)03]o.68-[PbTi03]o.32 (011) using surface-sensitive 
x-ray magnetic circular dichroism spectroscopy at the Co L 3,2 edges. Three distinct electric held 
driven remanent magnetization states can be set in the Co him at room temperature. Ab-initio 
density functional theory calculations unravel the relative contributions of both strain and charge to 
the observed magnetic anisotropy changes illustrating magnetoelectro-elastic coupling at artihcial 
multiferroic interfaces. 


Multiferroic systems, where two or more ferroic prop¬ 
erties (ferromagnetism, ferroelectricity or ferroelastic- 
ity) coexist, provide the opportunity to study coupling 
mechanisms between different order parameters [1]. The 
prospect of electric field control of magnetism with its 
potential use in technological applications ma has fo¬ 
cused attention on the subgroup of multiferroics that ex¬ 
hibits magnetoelectric (ME) coupling. Heterostructures 
consisting of cross-coupled ferromagnetic (FM) and fer¬ 
roelectric (FE) layers [5] are often referred to as artificial 
multiferroic composites. Due to their modular nature 
the number of available systems with potential multifer¬ 
roic properties is greatly increased compared to intrinsic 
multiferroic systems which proves to be an advantageous 
concept to achieve electric field control of magnetism at 
room temperature. 

The mechanisms involved in ME interface coupling of¬ 
ten result from new and interesting underlying physical 
phenomena. Strain-coupled systems [7HS1 make use of 
the piezoelectric properties of a ferroelectric system to 
control the magnetism in a ferromagnet through magne¬ 
tostriction. Furthermore, ferroelectric polarization rever¬ 
sal may change the overlap between atomic orbitals at the 
FM/FE interface [TUm^ or drive a magnetic reconstruc¬ 
tion at the interface. Charge-mediated ME cou¬ 

pling exploits the electric field effect [TM^ as well as the 
remanent electric polarization of FE components [23l [24] 
to modulate the charge carrier concentration in an adja¬ 
cent FM layer, where accumulation or depletion of spin- 
polarized electrons results in a change of the interface 
magnetization. Different length scales apply to the afore¬ 
mentioned mechanisms. While the influence of strain ex¬ 
tends to several tens of nanometers, the screening of sur¬ 
face charge takes place within the Thomas-Fermi screen¬ 
ing length (on the order of a few Angstroms in metals 


[IS]). The coexistence of strain and charge effects have 
seldom been reported [251(55] and so far been explained 
in a phenomenological framework. In this work, we dis¬ 
entangle strain and charge contributions to the magnetic 
response upon electrical switching using surface-sensitive 
x-ray magnetic circular dichroism (XMCD) at the Co 
L 3,2 edges and ab-initio density functional theory (DFT). 
The heterostructure consisting of a Co wedge (0-7 nm) 
grown on top of the ferroelectric [Pb(Mgi/ 3 Nb 2 / 3 ) 03 ]o .68 
-[PbTi03]o.32 (Oil) (from here on PMN-PT) allows for a 
thickness-dependent study. We find experimentally that 
it is possible to set three distinct remanent and reversible 
magnetization states through magnetoelectro-elastic con¬ 
trol at room temperature. DFT calculations for differ¬ 
ent strain and charge states reproduce the experimental 
behavior and unravel the different mechanistic contribu¬ 
tions. 

Relaxor FE PMN(- 3 _ 3 .)-PT 2 ; (Oil), with a composition 
of a; = 0.32 located in the morphotropic phase boundary 
region [25], (Atom Optics Co., LTD., Shanghai, China) is 
used as a substrate due to its strong piezoelectric proper¬ 
ties. Its crystal structure is monoclinic with lattice con¬ 
stants 0=4.02 A , 6=4.01 A and c=4.03 A [29]. A cobalt 
wedge with linearly increasing thickness from 0 — 7 nm is 
grown on PMN-PT (Oil) via thermal evaporation. X-ray 
diffraction showed that Co grows with fee (lll)-texture. 
A capping of 2 nm Cr prevents oxidation and a 30 nm Au 
film serves as bottom electrode. Figure (a) shows the 
sample design and measurement geometry. Depending 
on the electric field applied across PMN-PT (Oil), three 
distinct remanent FE polarization states can be set. The 
FE polarization is poled positively or negatively out of 
plane (OOP-1- or OOP-) by applying an electric field of 
±0.36 MV/m at the bottom electrode while the top elec¬ 
trode is connected to ground. When comparing OOP± 
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and OOP- poled FE no lattice parameter change in PMN- 
PT is expected and the Co top layer encounters identical 
strain conditions. However, FE polarization switching 
alters the interfacial charge that has to be screened by 
the adjacent cobalt layer through accumulation or de¬ 
pletion of electrons. Sweeping between opposite OOP 
FE polarization directions, PMN-PT (Oil) exhibits a re¬ 
manent in-plane (IP) poled state at the coercive elec¬ 
tric field (±0.14 MV/m). The switching from an OOP 
to an IP poled configuration and vice versa is accompa¬ 
nied by structural changes of the PMN-PT [301131] as 
indicated in Fig. (b) and (c) that act on the Co top 
layer. OOPoIP switching alters both the strained state 
of cobalt and the interfacial charge seen by the Co film. 
Note that both OOP poled states as well as the IP poled 
configuration are stable at remanence. The FE polar¬ 
ization of PMN-PT (Oil) at 298 K was measured to be 
2 • FVmn-pt =60 /iC/cm^. 

XMCD [32] measurements at the Co ± 3,2 edges 
were carried out at the X-Treme beamline |33| at the 
Swiss Light Source, Paul Scherrer Institut, Switzer¬ 
land and at beamline 6.3.1 [31] at the Advanced Light 
Source, Lawrence Berkeley National Laboratory, Califor¬ 
nia, USA. In XMCD, the absorption intensity difference 
between opposite light helicities is an element sensitive 
probe of magnetization along the photon propagation 
direction. Spectra were recorded at room temperature 
in grazing incidence geometry, where x-rays are inclined 
from the surface plane by an angle of 30°, measuring 
the projected magnetization along the (Oil) crystal di¬ 
rection of the PMN-PT. The external magnetic field was 
applied along the x-ray beam direction. Sum rules allow 
for a quantification of the Co spin and orbital magnetic 
moments and morb from analysis of the XMCD 

spectra [35] [36]. 

Co XMCD hysteresis loops along the (Oil) direction, 
taken in total electron yield (TEY) mode for oppositely 
OOP poled states as well as the IP poled configuration 
at distinct thicknesses of the wedge, highlight two differ¬ 
ent ME coupling mechanisms at play. Electrical switch¬ 
ing from an OOP poled to an IP poled state induces an 
anisotropy change with higher remanent magnetization 
as seen in Fig.[^(d) for a nominal Co thickness of 3.5 nm. 
The same behavior is observed probing a thicker part of 
the wedge at a nominal Co thickness of 6.3 nm in Pig. [^ 
(e). Additionally we observe a more subtle anisotropy 
change comparing hysteresis curves taken for oppositely 
OOP poled PE in Fig. (f). This anisotropy change is 
not observed in the thicker part of the wedge, as seen by 
the nearly identical hysteresis loops of Fig. (g). TEY is 
a surface-sensitive detection mode where the probability 
of electron escape from the Co/PMN-PT interface decays 
exponentially with increasing Co top layer thickness (the 
electron sampling depth for Co is about 2.5 nm [37]). 
Therefore, the observed difference in magnetic anisotropy 
in Pig. (f) and its absence in Fig. (g) hints that its 
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FIG. 1: (a) Measurement geometry. (b,c) Lattice parame¬ 
ter changes in OOP/IP poled PMN-PT, respectively. Green 
(red) arrows indicate compressive (tensile) strain in Go. (d-g) 
XMGD hysteresis curves probing the Co magnetization pro¬ 
jection along the (Oil) PMN-PT crystal direction for the three 
distinct FE polarization states. (d,e) Switching the FE po¬ 
larization from an OOP (blue curve) to an IP poled state 
(red curve) probing a nominal Co thickness of 3.5 nm (d) and 
6.3 nm (e), induces an anisotropy change with higher rema¬ 
nent magnetization, (f) For 3.5 nm Co thickness, OOP poled 
polarization directions exhibit also different anisotropies, (g) 
For 6.3 nm Co thickness the anisotropy change for oppositely 
OOP poled FE is now absent. 


origin lies at the interface between Co and PMN-PT. As 
pointed out above, this effect cannot be attributed to 
a piezoelectric-magnetostrictive coupling since the struc¬ 
ture of PMN-PT in the two states is equivalent. Hence, 
this anisotropy change due to the substrates’ opposite 
OOP polarities suggests a charge driven magnetoelectric 
coupling. The anisotropy change shown in Fig.[^(d) and 
Fig. 0 (e) at both the thinner and the thicker part of the 
wedge can be understood in terms of the magnetostric¬ 
tion of cobalt in response to the lattice parameter changes 
of PMN-PT [31]. Since strain is a ‘bulk’ effect, its influ¬ 
ence persists throughout the whole Co film thickness. For 
a quantitative analysis, a series of XMCD spectra was 
taken as a function of applied electric field on the thin 
part of the wedge at 3.5 nm Co thickness at magnetic 
remanence after saturation in 2 T in total fluorescence 
yield (TFY). Sum rule analysis was used to extract the 
magnetic moment rritot = Ws,eff + "m-orh projected along 
the (Oil) direction (for details, see ‘Supplemental Mate¬ 
rial’). The resulting dependence on the electric field is 
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FIG. 2: Total magnetic moment along the (Oil) direction 
at remanence obtained from sum rule analysis as a function 
of applied electric field. The grey curve is a guide for the eye 
to link subsequent measurements. The arrow indicates the 
sweeping direction. Black squares (open circles) indicate mea¬ 
surements coming from OOP- (OOP-f) poling. The dashed 
red branches highlight the dependence of rritot on the FE OOP 
substrate polarity when switching between OOP-f and OOP-. 


given in Fig.[^ where the gray curve links successive mea¬ 
surements. rritot is strongest at the coercive electric held, 
where the FE polarization is rotated in-plane. Compar¬ 
ing measurements of oppositely poled FE, OOP- poled 
PMN-PT results in a smaller Co rritot than OOP-1- poled 
PMN-PT. Here, rritot depends solely on the FE polariza¬ 
tion state that the PMN-PT has been set in, irrespective 
of an actively applied bias voltage. Note that in 2 T 
applied held no dependence of the saturation magneti¬ 
zation on the FE polarization can be observed. At 2 T 
held applied along the easy (100) direction, the effective 
spin moment ms,eff = 1.64±0.16 fj,B and orbital moment 
rHorb = 0.131 ± 0.002 compare well with literature 
values [Ml 1^ . 

The impact of the FE order of PMN-PT on the elec¬ 
tronic and atomic structure of a Co top layer is twofold. 
We observe a hysteretic behavior of remanent rritot for 
OOP-b and OOP- poled FE suggesting a charge-driven 
magnetoelectric coupling contribution due to accumula¬ 
tion and depletion of electrons at the FM/FE interface. 
The contribution of charge to the change in total mag¬ 
netic moment is highlighted by the dashed red branches 
in Fig. Deviations occur only at the coercive electric 
held, where strain dominates while no net surface charge 
should be present. As the total moment at 2 T does 
not appreciably change with FE polarization but there 
is a signihcant change to rritot at magnetic remanence, 
we attribute these changes in magnetization to changes 
in effective magnetic anisotropy energy (MAE) of the Co 
him. To investigate the separate inhuences of strain and 
screening charge on the MAE we perform hrst-principles 
DFT calculations of bulk fee cobalt with each perturba¬ 
tion applied separately. 

For Co hlms thicker than 1.5 nm [MMT] . the shape 
anisotropy dominates the MAE and dictates an isotropic 


in-plane magnetization. This isotropy within the him 
plane is subsequently lifted by other MAE contributions. 
The bulk magnetocrystalline anisotropy for fee Co fa¬ 
vors an easy axis along the [111]- and equivalent cubic 
directions. However, for a (111) him the strong shape 
anisotropy disfavors the low energy crystalline directions. 
Moreover, the volume magnetocrystalline anisotropy is 
isotropic within the (111) him plane and thus creates no 
anisotropy even if its magnitude is altered. 

Another contribution to the MAE is magnetoelastic¬ 
ity, which exhibits lower order terms of the directional 
magnetization expansion |39j that are coupled to strain 
tensor elements {cij). For cubic symmetry its energy con¬ 
tribution is 

Emag-el = -\- 2622^2 + £ 330 ^ 3 ) 

+2B2{e230t2cr3 + ei3aia3 + ei2aia2) , 

where Bi are the cubic magnetoelastic constants and a 
is the corresponding direction cosine of the magnetiza¬ 
tion. For the [lllj-oriented fee Co him we transform this 
expression [39| (see ‘Supplemental Material’) into hexag¬ 
onal coordinates to yield, for the him plane magnetiza¬ 
tion: 

Em,ag—el,hex{4’) — ~'^{Bi-\-2B2){£iqq — €q^j) sin (0) (1) 

where e[ are the strain elements in the hlm-plane labeled 
with respect to the PMN-PT substrate and (p is the an¬ 
gle of the in-plane magnetization relative to the [100] di¬ 
rection. The magnetoelasticity creates an easy in-plane 
direction which is determined by an ‘effective’ magnetoe¬ 
lastic constant Heff = Bi + 2 B 2 ■ 

By performing total energy calculations for a set of 
strained fcc-cobalt unit cells (see ‘Supplemental Ma¬ 
terial’) we compute Bi and B 2 using DFT. We hnd 
both Bi=—8.7 MJm“^ and 52=7.2 MJm“^ in reason¬ 
able agreement with experimental and theoretical liter¬ 
ature values [3911421 [43]. Moreover, the combination of 
these values gives a positive effective magnetoelastic con¬ 
stant, 5eff. Consequently, we predict that a net strain 
(e'loo — Eg^j) > 0 creates an easy axis along the [Oil] di¬ 
rection, whereas (eioo ~ ^oiT^ ^ produce an easy 

axis parallel to [100]. In PMN-PT, OOP—)-IP poling is 
accompanied by a strong positive e'^gg transferred to the 
Co him [M] resulting in a positive net strain. Hence, 
our theoretical hnding is in agreement with the experi¬ 
mentally observed anisotropy change along [011] upon IP 
poling. 

For both the OOP-1- and the OOP- poled state, the Co 
him encounters a net strain (eioo ~ ^oiT^ ® 
perimentally observed magnetization shows a preferred 
orientation close to the [100] axis in agreement with our 
prediction. However, in the experiment there is a 15% 
higher magnetization projection along the [Oil] axis for 
the OOP-k state than for the OOP- state. Since the 
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structure of PMN-PT in the two states is equivalent, the 
difference has to be attributed to a contribution stem¬ 
ming from the FE polarization direction. 

For example, the presence of interface charge aint may 
necessitate screening by the valence electrons of the adja¬ 
cent Co film. With 2 • Ppmn-pt=60 /rC/cm^, the amount 
of interface charge doping for fee (111) Co can be esti¬ 
mated to be CTini(O) = ±0.102 e“/unit cell area. This 
charging will be largest at the interface and then de¬ 
cay exponentially corresponding to the Thomas-Fermi 
screening as aint{z) = where z measures 

the distance from the interface and Aco is the Thomas- 
Fermi-screening length of Co (Aco =0.15nm [H]). 

Next we examine the impact of this interface charge on 
the magnetoelastic constants (i?i((T), i? 2 (o')), as shown 
in Fig. ia), by repeating our computations with a var¬ 
ied total e~ count within the DFT calculations. We find 
(Fig.j^a)) a strong variation of Bi with charging whereas 
B 2 remains nearly unchanged. Moreover, the different 
behavior of Bi and B 2 as a function of charging leads to 
a sign change of B^g around cr =4:^CIcm‘^, as depicted 
by a dashed green line in Fig. |^a). For the same neg¬ 
ative net strain corresponding to OOP poled PMN-PT, 
the OOP± and OOP- cases have different alignments of 
the magnetic easy axis at the interface, as sketched in 
Fig. (b,d). For the OOP- case, the accumulation of 
positive charges at the interface {a < 0, BeS > 0) creates 
an easy axis along the [100] direction. On the other hand, 
in the OOP± state the accumulation of negative charges 
(cr > 0) reverses the sign of Seff and thus favors the or¬ 
thogonal [oil] direction as easy magnetization direction. 
Consequently, we expect that switching of the electric 
polarization in combination with an alternation of the 
magnetoelastic constants by interface charging leads to a 
90°change of the preferred magnetization direction. This 
magnetoelectro-elastic effect will be constrained to the 
interface region, where enough charge accumulation is 
present. 

The experimentally observed higher remanent magne¬ 
tization along the [011] direction for the OOP± state 
compared to the OOP- state in Fig. [^ and reflects 
our calculated magnetoelectro-elastic effect. However, 
the detected signal contains contributions of both, the 
strain that extends throughout the entire film, as well as 
the charge — an interface effect. Consequently, the ex¬ 
ponential decay of the charge screening away from the 
Co/PMN-PT interface suppresses the measurement of 
the predicted effect in thicker films when using surface- 
sensitive TEY detection mode (Fig. [^g)). 

In conclusion, we investigated by a combined experi¬ 
mental and theoretical effort the magnetic properties of 
the artificial multiferroic Co/PMN-PT interface. From 
our XMCD measurements we found that the magnetic 
anisotropy behavior of the Co film depends on the three 
distinct polarization states (IP, OOP(±,-)) the PMN-PT 
can be set in. For thin film thicknesses in which inter- 




FIG. 3: (a) Variation of magnetoelastic constants as a func¬ 
tion of charging of the unit cell cr. (b,c,d) Sketch of the three 
switching states (OOP-, IP, OOP-f) of the PMN-PT/Co in¬ 
terface. The black/red arrows show the direction of elec¬ 
tric/magnetic polarization and the blue arrows show the net 
strain, (b) illustrates the OOP- case where strain and elec¬ 
tron charge lead to a preference of the [100] direction as the 
easy axis, (c) depicts the highly strained IP state with no 
interface charging and an easy axis along [011]. In (d) the 
combination of strain and positive interface charge creates an 
easy axis along [Oil] at the interface which decays and turns 
towards [100] away from the interface. 


face effects dominate we find a significant difference be¬ 
tween all three states, whereas for thick film thicknesses 
the difference between the OOP states vanishes. Our 
theoretical investigation illustrates that the changes in 
anisotropy are due to a combination of magnetoelastic¬ 
ity and interface charging opening up the possibility for 
enhanced magnetoelectric coupling. Finally, we suggest 
that the found modulation of magnetic anisotropy by the 
magnetoelectro-elastic effect may allow to create a mag¬ 
netic anisotropy gradient in thin films. If the gradient is 
strong enough, it could give rise to a spiral state in the 
thin film, which could be controlled by the ferroelectric 
substrate. 
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